Many-Body Coherent Destruction of Tunneling 
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A new route to coherent destruction of tunneling is established by considering a monochromatic 
fast modulation of the self-interaction strength of a many-boson system. The modulation can be 
tuned such that only an arbitrarily, a priori prescribed number of particles are allowed to tunnel. 
The associated tunneling dynamics is sensitive to the odd or even nature of the number of bosons. 
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The phenomenon of coherent destruction of tunneling 
(CDT) [H, 0] by a driving field has been one seminal re- 
sult in studies of quantum dynamics control. Direct ob- 
servation of CDT was recently achieved in several exper- 
iments [3, 01 , one of which [j] involving non- interacting 
cold atoms in a double-well potential. CDT in interact- 
ing many-body systems has also attracted considerable 
interest, with previous work focusing on driving fields on 
resonance with the interaction energy 0, 0] . The Mott- 
superfluid transition in ultracold systems via a mecha- 
nism similar to single-particle CDT has also been ob- 
served experimentally [7|. 

Traditionally, the driving field in CDT studies is to 
directly modulate the bare single-particle levels of an 
undriven system. For example, in a two-level theory of 
CDT, the driving field is to modulate the energy diff'er- 
ence between two bare levels. By contrast, in this work 
we expose a new route to CDT by taking advantage of 
the particle-particle interaction in a many-body system. 
Specifically, we consider the CDT in a two-mode Bose- 
Hubbard model that describes a two-mode Bose-Einstein 
condensate (BEC). We show, both analytically and com- 
putationally, that a monochromatic off-resonance driving 
of the self-interaction strength of the BEC can induce 
different types of CDT, without a direct modulation of 
the mode-energy bias. Interestingly, this makes it pos- 
sible to precisely control, at least in principle, the num- 
ber of bosons allowed to tunnel. Another remarkable 
prediction is the sensitivity of the full-quantum dynam- 
ics to the even or odd nature of the number of bosons. 
Note that in other contexts such as matter-wave solitons, 
intriguing effects of a periodic modulation of the self- 
interaction strength of a BEC (the so-called "Feshbach- 
resonance management") have been discovered Q, but 
on the mean-field level only. 

Consider then the following Bose-Hubbard Hamilto- 
nian for a two- mode BEC, 

H = vh{a\ar + alai)/2 + g{t)'h{a\ai - a3'.a,.)^/4, (1) 

where r and I are mode indices, au and a\. [k — r, I) 
are the bosonic annihilation and creation operators, v 



describes the constant tunneling rate between the two 
modes, and g{t) is the interaction strength between same- 
mode bosons. We use the unit of v to appropriately scale 
all the parameters such that v, g{t), and t all become 
dimensionless variables. The total number of bosons 
N — ajai + alttr is a conserved quantity and the dimen- 
sion of the Hilbert space is iV + 1. Using the Schwinger 
representation of angular momentum operators, namely, 
Ja; = {ajar + alai)/2, Jy = {alai — a\ar)/{2i), and 
Jz = {ajai — alar)/2, eq. ^ reduces to 



H{t) = vTiJ^ + g{t)njl 



(2) 



The Hilbert space is expanded by the eigenstates of Jz, 
denoted |m), with Jz\'m) = to|to). The mode population 
difference is given by the expectation value of 2Jz. For 
later use we also define Jj^ = Jx + iJy Note that there is 
no energy bias between the two modes, and that the time 
dependence of the Hamiltonian arises from g{t), which is 
assumed to be Isl 



9{t) =50 + 51 cos(wt). 



(3) 



Under appropriate conditions H{t) describes a BEC dis- 
tributed in the two wells of a double- well potential 01, 
in the ground band and the first-excited band associated 
with an accelerating optical lattice or in two hyper- 
fine levels. For convenience, below we focus on the first 
context, which can be best realized by optical superlat- 
tices [lOj. Hence r{l) denotes the right (left) well. Our 
central idea is to use an off-resonance oscillation in g{t) to 
switch off the tunneling between the left and right wells. 
The Floquet operator associated with H{t) is given 

by F = r [exp(-i J^''^'^ H{t')/h dt')], where T is the 
time-ordering operator. Its eigenstates are the Floquet 
states, with eigenvalues exp(— ie27r/?iw), where e is the 
quasi-energy. Because H{t) apparently possesses a par- 
ity symmetry, namely, it is invariant upon an exchange 
of the indices I and r, the Floquet states can be chosen 
as either positive-parity or negative-parity states. If two 
opposite-parity Floquet states cross or touch, then simi- 
lar to the single-particle CDT mechanism, their superpo- 
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FIG. 1: (Color online) Quasienergy spectrum versus 
gi/io. Vertical dashed lines indicate the location (i.e. 
gi/uj ~ 0.267, 0.344, 0.481) of level degeneracies. Dotted 
(solid) lines are for states with negative (positive) parity. 
iV = 10, go = 0, and v = 1. In (a) a; = 40 and in (b) 
uo — 2tt. Spectral details around the square are shown in in- 
set in (b), with two level-crossings. Here and in other figures 
all variables are dimensionless. 



sition, which is still a Floquet state, breaks the left-right 
symmetry and hence CDT occurs 
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Figure 1 depicts that the degeneracy between opposite- 
parity states can easily occur, thus suggesting that CDT 
is possible via solely a monochromatic modulation in the 
self-interaction strength. Next we examine some detailed 
features presented in Fig. 1, for go = 0, N = 10, and 
ti = 1 as an example. In Fig. l{a.), uj — iO ^ v ^ 1, 
and three values of gi/to for which level degeneracies oc- 
cur are explicitly marked by vertical dashed lines. At 
point I, the level crossing involves two positive-parity 
states and one negative-parity state. At point II, two 
pairs of opposite-parity states become degenerate simul- 
taneously. At point III, three states in the middle of 
the Floquet spectrum cross, and two additional pairs 
of opposite-parity Floquet states also touch each other. 
Similar spectral patterns can be found in other regimes 
of gi/oj. In Fig. 1(b), = 27r is quite comparable to 
V = 1. In this intermediate- frequency case the spectral 
patterns become more complicated. Nevertheless, as in- 
dicated by those vertical dashed lines in Fig. 1(b), when 
gi increases, the spectral pattern in Fig. 1(b) becomes 
analogous to those seen in Fig. 1(a). 

We develop below a theory for the Floquet spectrum 



in the high-frequency regime, with its validity condition 
elaborated later. In that regime, the full Floquet the- 
ory can be expanded to the first order of 1/w [3] and a 
static effective Hamiltonian HcS for the driven quantum 
dynamics can be obtained by averaging out the driving- 
field effects Ol. Explicitly, 



2^ 



2tt/u 



where A{t) — J^giCos{ujt) dt — (gi/cj) sin(a;t). 
Using the identity in the SU(2) algebra, i.e., 
e^^W-^J J^e-^^W-^' = (J+/2)e*^W(2J.+i)_^p_c_^^ijgj.gp,,^ 

means complex conjugate of the preceding term, and sub- 
stituting this into Eq. (jH) to perform the integral, we 
obtain the effective Hamiltonian, 

Hcff = gohJl + {vhJ+/2)Jo [.gi(2J, + l)/uj] + c.c, (5) 

where Jq{x) is the ordinary Bessel function of order zero. 
Equation ([5]) indicates that the net effect of a fast mod- 
ulation in g{t) is the rescaling factor Jo[5i(2^z + 
which depends on J^, i.e., the population difference be- 
tween the two wells. Further, a nonzero go leads to the 
go^Jz term in HcS- It is well known that this term can 
induce population localization via a self-trapping mech- 
anism. To isolate population localization due to possible 
CDT phenomena from that due to self-trapping, we will 
not consider cases with nonzero until much later. 

In the eigen-representation of Jz, Hcs is a tri-diagonal 
matrix. Further, if (m— Iji/cfflm) = 0, then {m\Hcs\m — 
1) = and we must also have (1 — mjiJcffl — = 
(— m|i7cff|l ^ Ti) = due to symmetry consideration. 
These four zero matrix elements divide the tri-diagonal 
matrix of Hcs into three uncoupled subspaces, i.e., 



HcS — 



hi 
h, 
hr 



(6) 



where hi represents a sub-matrix of Hcs in the subspace 
spanned by states \N/2), \N/2 — 1), • • • , |m) (assuming 
m > 0) , hr represents a sub-matrix of Hcs in the sub- 
space spanned by states | — N/2) , | — N/2 + 1) , • • • , \ — ni), 
and hi represents the third block matrix involving other 
remaining basis states. Therefore, if (to — l|i7cff|TO) — 0, 
then (i) the transition between \m) (| — to)) and all other 
basis states |to') (| — to')) with m' < m will not occur; (ii) 
the Floquet states must display degeneracy because hr is 
identical with hi due to the parity symmetry of H{t); and 
(iii) the dimension of hr or hi, namely, {N/2 — m -f 1), 
also gives the expected number of degenerate pairs. 

Suppose there are N — i {i) particles in the left (right) 
well. Without loss of generality we assume i < N/2. 
The associated quantum state is given by \m) — \N/2 — 
i). Using Eq. (O, one finds that the condition (m — 
l\Hcfi\m) = is equivalent to 



Jo[gi{N ~{2i + l))/u]=0. 



(7) 
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FIG. 2; (Color online) (a) Time dependence of (S) for go — 
and A'' = 10, with = 40 in (a) and = 27r in (b). From 
top of bottom, values of gi/uJ are associated with the three 
vertical lines (from left to right) in Fig. [ija) for panel (a) 
and in Fig. [ijb) for panel (b). For comparison, a case with 
Qi/io — 10.5/40= 0.2625 is also shown in (a). Inset of (b) is 
for the first level crossing shown in the inset of Fig. [Hb) . 



As such, if gi{N — (2« + l))/ti; is tuned to become one 
root of j7o(a;), then the tunneling of one more particle 
from the left well to the right well (hence m ^ m — 1) 
becomes prohibited if the left well has already released i 
particles to the right well. 

We now compare the theoretical result of Eq. ^ with 
our computationally precise results in Fig. 1. For the 
cases in Fig. 1(a), theoretically the spectral degeneracy 
is expected to occur when 51 [10 — {2i + l)]/oj becomes a 
root of Ja{x). In particular, for the first root of Jo{x) 
at X ^ 2.405 and for i — 0,1,2, the predicted degener- 
acy is at gi/uj ~ 0.267,0.344 and 0.481, with the pair 
number of level degeneracies in each case given by the 
dimension of hr, i.e., i + 1. This is in perfect agreement 
with the three marked degeneracy points shown in Fig. 
1(a). It is now also possible to explain why the level 
crossings in the middle of the spectrum in Fig. 1(a) can 
involve three states. This is because: (i) if N is even 
and go — 0, i?eff always has a zero eigenvalue due to 
its tri-diagoal structure; and (ii) when the dimension of 
and hi is odd, they can present two additional zero 
eigenvalues. Interestingly, for sufficiently large t/i, such 
agreement between theory and numerics may persist for 
the intermediate-frequency case in Fig. 1(b). For exam- 
ple, the degeneracies marked by the three vertical lines in 
Fig. 1(b) occur at g^juj - 2.007, 2.133, and 2.986. These 



values, when multiplied by (N — {2i + 1)) for i = 0, 1, 2, 
respectively, are the 6-th or 5-th root of j7o(x). Note 
however, that the subtle crossing behavior depicted in 
the inset of Fig. 1 (b) is beyond our theory. 

With a normalized population imbalance (S) = 
2{Jz)/N, Fig. IHa) shows the numerically exact popu- 
lation dynamics associated with the three marked points 
in Fig. 1(a). The initial state is that all particles are in 
the left well. In the first case for gi/uj ~ 0.267, (5*) is seen 
to stay at almost unity and hence in essence the tunneling 
between the two wells is completely suppressed. In the 
second case for gi/uj ~ 0.344, our theory predicts that 
the tunneling suppression occurs only when (S) becomes 
0.8. As seen in Fig. 2(a), (5*) indeed oscillates between 
0.8 and 1.0. Similarly, in the third case, (S) oscillates 
between 1.0 and ~ 0.6, confirming our theory that the 
tunneling stops if two particles are already released to 
the right well. Excellent agreement is obtained at other 
level degeneracy points. These features signify one key 
aspect of our many-body CDT: it depends sensitively on 
the number of particles that have already tunneled. With 
the same initial condition, in Fig. 2(b) we also show the 
three intermediate-frequency cases marked earlier by the 
vertical lines in Fig. 1(b). The associated population dy- 
namics still agrees with our theory. Because the predicted 
CDT points are independent of the actual tunneling rate 
V, we found that even if an oscillation in v is considered 
(which can be induced by the modulation in g{t)), anal- 
ogous results can be obtained. For the case shown in the 
inset of Fig. 1(b), which is beyond our high-frequency 
theory, almost complete population delocalization is ob- 
served in the inset of Fig. 2(b). 

Figure 3(a)-(c) schematically illustrates the three rep- 
resentative CDT cases studied in Fig. 1(a) and Fig. 2(a). 
In Fig. 3(a), no particle is allowed to tunnel; and in 
Fig. 3(b) and 3(c) one or two particles have tunneled 
and then CDT occurs. Consistent with this picture. Fig. 
3(d) depicts the numerically time- averaged (5*), denoted 
((5)), as a function of gi/u), for = 10 and uj = 40. It 
is seen that as gi/to is scanned at a rather low resolu- 
tion, the value of {{S)) is either zero or close to "magic" 
nonzero numbers 1.0, ~ 0.9, and ~ 0.8, • • • ). Figure 
3(a)-3(c) also provoke us to re-interpret our theoretical 
finding. In particular, for go = 0, the energy difference 
between the two configurations in Fig. 3(a) and Fig. 3(b) 
(Fig. 3(b) and Fig. 3(c)) is gihcos{ujt)[N - (2i + 1)] 
with i = (i = 1). As such, even though there is 
no direct modulation of the energy bias between the 
two wells, the oscillation in g{t) still causes a modu- 
lation of the effective bias between different configura- 
tions. With this interpretation we are able to re-derive 
Eq. I?]) by analogy to standard single-particle CDT 
theory. This analogy also makes clear that the precise 
condition for our high-frequency approximation should 
be w ^ max(w,Y^) 0], where v = u ■\/(]V^^7)(i"TT) 
is the coupling strength between states \N/2 — i) and 
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FIG. 3: (Color online) (a)-(c): Schematic picture of particle 
number-dependent CDT. In case (a) no particle tunnels, and 
in cases (b) and (c) only one or two particles tunnels and 
then CDT occurs, (d): Long-time average of ((S)) vs gxju) 
for tj = 40 and go = 0, with gijio scanned at a rather low 
resolution, (e): Same as in (d), but with gi/cJ scanned in 
small steps around the regime of 0.343, with go/w = 1/144 
(squares), go/uo = 1/288 (triangles), and po = (circles). The 
total time used for averaging is 20000 in dimensionless units. 



\N/2 - i - 1), and e = \gi[N - {2i + 1)]| is the amplitude 
of the effective bias. For i ^ N, this condition becomes 
V <^ max ( , y^giuj) . This well explains our early ob- 
servation from Fig. 1(b) that as 171 increases, the spectral 
pattern in intermediate-frequency cases starts to resem- 
ble those in Fig. 1(a) and becomes more perspicuous 
with our theory. 

Let us now turn to cases with nonzero go- If the sys- 
tem is close to a CDT point, Jo [giC^Jz + 1)/'^] is small, 
the go term in Hcs will dominate and hence the associ- 
ated self-trapping effect may induce a strong population 
imbalance on a very long time scale. Taking one small 
window of 51/1^ in Fig. 3(d) as an example, we com- 
pare 50 7^ with go — cases in Fig. 3(e). Clearly, 
as go increases, the width of the {{S)) profile in Fig. 
3(e) increases significantly. This interplay between self- 
trapping and CDT is analogous to that in a two-mode 
optical waveguide system where the mode bias is period- 
ically modulated [l3|. The peak value of {{S)) is also seen 
to change with go ■ We conclude that on one hand a small 
nonzero go is beneficial to experiments because it reduces 
the sensitivity of {{S)) to the exact values of gi/oj; on the 
other hand, however, the predicted particle-number de- 
pendent CDT effect may be buried by self-trapping if go 
is too large. 

Our findings have potential applications in probing and 
exploring genuine quantum coherence in BEC. In particu- 
lar, the dynamics under a pre-established CDT condition 
may be dramatically changed upon adding bosons to the 
system. As an example we consider the CDT point I in 



Fig. 1(a), where Jo [giiN - l)/w] for 10. If we 
let ^ -I- 2 by adding two particles to the left well, 
then because we still have jTo + 2) — 3)/w] = 0, 

Eq. ([7]) suggests that the CDT will be re-established af- 
ter one particle is tunneled to the right well. However, 
if we let A^ — > A^ -I- 1 by adding only one particle to 
the left well, then because (A^ — 1) cannot be written as 
(iV + 1) - (2i + 1) for any i, Jo [giUN + 1) ^ {2i + l))/uj] 
is in general nonzero for fixed gi/ui, and as a result all the 
particles start to tunnel back and forth between the two 
wells. The minor difference between adding an even and 
adding an odd number of bosons is thus greatly amplified 
by CDT, a prediction also confirmed by our numerical ex- 
periments. Similar behavior is obtained if more particles 
are added to the system. This odd-even sensitivity to the 
particle number is absent in any mean-field theory of a 
BEC, providing a possible means for accurate counting 
or efficient filtering of the number of bosons. Certainly, 
as implied by the results in Fig. 3(e), this is possible 
only if go is sufficiently small such that the {{S)) profiles 
associated with different particle numbers do not overlap. 
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